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The Riccati formulation of the quasiclassical theory of nonequilibrium superconductors is devel- 
oped for spin-dependent scattering near magnetic interfaces. We derive boundary conditions for the 
Riccati distribution functions at a spin-active interface. The boundary conditions are formulated in 
terms of an interface S-matrix describing the reflection and transmission of the normal-state conduc- 
tion electrons by the interface. The S-matrix describes the effects of spin-filtering and spin-mixing 
(spin-rotation) by a ferromagnetic interface. The boundary conditions for the Riccati equations are 
applicable to a wide range of nonequilibrium transport properties of hybrid systems of supercon- 
ducting and magnetic materials. As an application we calculate the spin and charge conductance 
of a normal metal-ferromagnet-superconductor (NFS) point contact; the spin-mixing angle that 
parameterizes the S-matrix is determined from experimental measurements of the peak in the sub- 
gap differential conductance of the NFS point contact. We also use the new boundary conditions to 
derive the effects of spin-mixing on the phase-modulated thermal conductance of a superconducting- 
ferromagnetic-superconducting (SFS) point contact. For high-transparency (metallic ferromagnet) 
"tt" junctions, the phase modulation of the thermal conductance is dramatically different from that 
of non-magnetic, "0" junctions. For low-transparency (insulating ferromagnet) SFS tunnel junctions 
with weak spin mixing resonant transmission of quasiparticles with energies just above the gap edge 
leads to an increase of the thermal conductance, compared to the normal-state conductance at Tc, 
over a broad temperature range when the superconducting phase bias is > tt/2. 



I. INTRODUCTION 

Spin dependent transport in hybrid systems composed 
of superconductors and spin-active materials such as fer- 
romagnets has attracted a lot of attention because of 
the possibility of generating coherent spin transport for 
spintronic devices.^ The spin polarization of a ferromag- 
netic material, one of the key parameters in the devel- 
opment of spintronic devices, is usually measured either 
by spin-dependent tunnelling techniques,^ or by point- 
contact Andreev reflection spectroscopy.^ Both methods 
infer information about the spin polarization from the 
conductance data of superconductor-ferromagnet (SF) 
junctions. When sandwiched between two s-wave su- 
perconducting leads (an SFS junction), a ferromagnetic 
layer can produce a "tt junction" , i.e. a ground state 
of an SFS junction in which there is a tt phase differ- 
ence between the two superconductors.^'^ The tt state 
has been observed in SFS junctions with metallic ferro- 
magnetic layers. It is predicted theoretically that an 
insulating or semiconducting ferromagnetic layer can also 
produce a tt junction. ^'^'^ More complicated tt junctions 
in which the Josephson coupling is provided by an inho- 
mogeneous magnetization, e.g. a ferromagnet-insulator- 
ferromagnet (FIF) tri-layer, have also been investigated 
theoretically. ^^^^^ 

The study of SFS junctions is fuelled in part by the 
proposal that tt junctions can be used to construct a 
non-dissipative superconducting phase qubit.^^ Most the- 
oretical investigations of SFS junctions are restricted to 
equilibrium properties, however the performance of a tt 
junction as a qubit depends sensitively on the suppres- 
sion of dissipative dynamics under nonequilibrium condi- 
tions. Recently, the nonequilibrium transport properties 
of Josephson junctions with spin active interfaces have 



begun to be explored theoretically.^^ 

A powerful formulism for calculating the nonequi- 
librium properties of superconducting heterostruc- 
tures is provided by the quasiclassical theory of 
superconductivity. Traditionally one obtains the 
quasiclassical Green's functions by solving the transport 
equations subject to boundary conditions at surfaces or 
interfaces. 

A multi-step approach to the boundary value problem 
at a surface or interface based on an interface transition 
matrix has been used by several authors. ^^■^^"^"'^ This 
method requires one to calculate an auxiliary Green's 
function for an impenetrable surface. The auxiliary 
Green's function is then used as an input to a T-matrix 
equation from which one constructs the quasiclassical 
Green's functions at the interface. This method can 
be applied to a broad class of interface models, and is 
suitable for numerical computations, ^^"^^ but it requires 
the computation of intermediate, unphysical auxiliary 
Green's functions. 

Boundary conditions which are expressed only in terms 
of the physical quasiclassical propagators and interface 
reflection and transmission amplitudes have been de- 
rived from microscopic scattering theory by Zaitsev^^ and 
Kieselmann^^ for nonmagnetic interfaces, and for spin ac- 
tive interfaces by Millis, Rainer and one of the authors. 
These boundary conditions are formulated as a set of 
third order equations in terms of the matrix Green's func- 
tions at the boundary, connected via an interface scatter- 
ing matrix (S-matrix) for normal metal electrodes. Al- 
though auxiliary propagators are not present, the non- 
linear boundary conditions are non-intuitive, difficult to 
solve and contain unphysical solutions which must be dis- 
carded. 

Recently, a more intuitive and computationally effi- 
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cicnt form of the quasiclassical boundary conditions was 
obtained for non-magnetic interfaces by Eschrig.^^ This 
formulation starts from the boundary condition of Za- 
itsev and Kieselmann, and is obtained by parameteriz- 
ing the quasiclassical Green's functions in terms of Ric- 
cati amplitudes. ^^"'^^ By formulating the quasiclassical 
theory in terms of the Riccati amplitudes, not only arc 
the transport equations easier to solve numerically, but 
the boundary conditions become linear and free of spu- 
rious solutions. Eschrig's formulation of the bound- 
ary condition amounts to finding physical solution to 
the Zaitsev-Kieselmann nonlinear boundary condition. 
For spin-active interfaces Fogelstrom obtained bound- 
ary conditions for the retarded and advanced coherence 
functions.^ However, a complete set of boundary condi- 
tions for non-equilibrium transport with spin-active in- 
terfaces was lacking. 

In this paper we derive the boundary condition for 
the quasiclassical Riccati amplitudes, both the coher- 
ence functions and distribution functions, for spin ac- 
tive interfaces and apply the new boundary conditions to 
study the nonequilibrium transport properties of clean 
superconductor-ferromagnet hybrid systems. The paper 
is organized as follows. The complete set of boundary 
conditions for the Riccati amplitudes at spin active in- 
terfaces is presented in section II, with technical steps 
of the derivation described in an appendix. In section 
III, the S-matrix for scattering by two models for spin- 
active interfaces, a ferromagnetic-insulating layer and a 
ferromagnetic- metallic layer, are derived and discussed in 
terms of the effects of spin Bltering and spin mixing. Ap- 
plications of the theory to the conductance of the normal 
mctal-ferromagnct-superconductor junction is analyzed 
in section IV. In section V the influence of spin mixing 
on the phase sensitive heat transport in SFS point con- 
tact is discussed in detail. 



II. THE BOUNDARY CONDITIONS FOR 
RICCATI AMPLITUDES 

In the Riccati formulation of the quasiclassical the- 
ory of nonequilibrium superconductivity, the quasipar- 
ticle excitation spectrum is determined from coherence 
functions, j^'^ and ^^M, which measure the relative 
amplitudes for normal-state quasi-particle and quasi-hole 
excitations; the occupation probability of theses states is 
described by distribution functions, and .^'^^^^ For 
brevity we refer to both types of functions as Riccati am- 
plitudes, or Riccati functions, since all obey Riccati-type 
transport equations, defined on classical trajectories in 
phase space (p,R). Thus, in general the Riccati ampli- 
tudes arc functions of space, R, time, t, the direction of 
the Fermi momentum, pf (or Fermi velocity vy) and the 
excitation energy, e. 

The Riccati amplitudes depend on spin, and in gen- 
eral are described by 2 x 2 density matrices in spin space 
whose eigenvalues determine the local coherence and dis- 



tribution functions for two possible spin states. The co- 
herence amplitudes obey Riccati-type equations; for ex- 
ample, 

ivy V7^^ = -2e7«+7«A7« + SV-7^S^-A^. (1) 

The distribution function, x^ , obeys a Riccati-type 
transport equation, 

z(at + v/-V)x^ = (7^A^ + E^).T^ 
+ :e^(A'47'4_s^) 

- 7«S^7^-S^ 

+ A^7^+7^A^, (2) 

where E'^ and A^, /i = R,A,K, are the diagonal and 
off-diagonal self-energies, respectively. Wc follow the no- 
tation in Ref. 27 throughout the paper. Particle«-+hole 
conjugation, denoted by ~ is defined by the operation 
q{pf,e) = q*{—pf,—e). The product of two fimctions 
of energy and time is defined by the non-commutative 
convolution, 

AB = Ao B{e,t) = e''-^'^^*''^^^^^^/'^Aie,t)Bie,t) . (3) 

Neither the operator, o, nor the arguments, (e,t), are 
shown explicitly unless required. 

Once the Riccati equations are solved subject to ap- 
propriate boundary conditions, the quasiclassical Green's 
functions can be constructed from the Riccati ampli- 
tudes. Physical observables such as the charge or heat 
current can then be calculated. This procedure is dis- 
cussed extensively by several authors, c.f. Refs. 27,30. 

At an interface or surface the local electronic poten- 
tial changes on an atomic length and energy scales. Such 
strong, short-range potentials are treated within the qua- 
siclassical theory as boundary conditions for the quasi- 
classical Green's functions, or equivalently the Riccati 
amplitudes. Such an interface can be described by a 
scattering matrix, S, for normal-state electrons and holes 
with excitation energies near the Fermi surface. Here 
we confine our discussion to specular interfaces, in which 
case the momentum of an excitation parallel to the in- 
terface, p||, is conserved. The interface S matrix is 
then described by a unitary matrix in the combined spin, 
particle-hole and direction spaces. Thus, 

SiPf) = ( I" t ) , (4) 

where the index 1 (2) refers to the left (right) side of the 
interface. Each element of this matrix, Sij, is a diagonal 
Nambu matrix in particle-hole space, 

in which Sij and S_ij are matrices in spin space, which 
are related by particle-hole conjugation, ^^^^(pii) = 
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FIG. 1: The Riccati amplitudes corresponding to the scatter- 
ing trajectories for a partially transmitting specular interface. 
The interface normal is z. The trajectories for the scattering 
states are labelled li, lo, 2i and 2o. 

P|[)]*' , where [..J'"^ is the matrix transpose in spin 
space. ° 

The Riccati amplitudes for the set of scattering 
trajectories labelled li, lo, 2i and 2o, shown in Fig. 
1, are classified into two groups. The quantities 
{if'^i Ij^'^ J J ^f}^ denoted by lower case symbols, are 
obtained by integrating the Riccati equations along the 
four trajectories from the bulk toward the interface. The 
quantities {T^'^, T^'"^, , denoted by upper case 

symbols, are obtained by starting at the interface and in- 
tegrating along the trajectory into the bulk. The bound- 
ary conditions listed below connect the unknown upper 
case amplitudes at the interface with the known lower 
case amplitudes via the interface scattering matrix. The 
derivation of these boundary conditions is outlined in the 
appendix. For example, the boundary conditions for tra- 
jectory lo are, 

rf = r^nMi + th2"s\,, (6) 

= ^ii7M+^i272^^il: (7) 

X^< - r^'r^P^ ^f^T'^fA -„R^K~A ,o\ 

where we have introduced effective reflection (r), trans- 
mission (t) and branch-conversion transmission (a) am- 
plitudes, 

r^i = +[i0?i)~'sl~iPg)-'sl,r\p^,)-\ (9) 

t?i = -[w^i)~'si~ {f3^,r'si,]-Hp^2r\ (10) 

= +ip2i)-'[Siiif3^i)-' ^ S,,if3^,)-T\ (11) 

it = -ip22)-'[SiiiPA)-' - S,,if3^,)-T\ (12) 

af, = (rf^ii-5n7f)(/3f2)-S (13) 

at = 0t2)-\S},,rt~l^Sl,). (14) 



The auxiliary quantities, P.^j , are defined as 

= 4 - 7f47f ; i3^,=S,^^ 7f %7f , (15) 

(3^ = s,,- its^.if; ^ = 3].^- its]at ■ (16) 

Similarly, for trajectory \i we have 

ff = f«7f5n+tf,72''52i, (17) 

rf = SWrt + Sll^tt, (18) 

= r'^^x^rt + t^ii^tt - af,xf af, , (19) 

where 

rli = +[(4'l)-'^ii-(/32l)-'^2i]-'(/32l)-', (20) 

= -{[fi^2)-^S^^^{M2r^S^^\-\M\)-\ (21) 

rt = +(/32l)"M5li(41)"^-4i(41)-']-S (22) 

tt = -(^2l)-MMi(/32l)-'-4i(/32^2)-']-S (23) 

~< - (ff^Ii-^L7f)(/3f2)-\ (24) 

4r = (/3i'2)-'(5iirf-7i^^n). (25) 

The boundary conditions for trajectories li and 2o are 
given by interchanging indices 1 ^ 2 in Eqs. (6)-(25). 
The derivation of Eqs. (6)- (25) is described in the ap- 
pendix. 

Note that there is more than one physically equivalent 
representation of the boundary condition for any of the 
coherence functions. For example, it is straightforward to 
show that an alternative form of the boundary condition 
in Eq. 6 for F^^ is given by 

rf = ^ii7iV. + 5i272''if. (26) 
rl = +0^2)-'[Sn0ur'-S,M2)-T' (27) 
t?r = -i^^2)-'[SuiP52)-'~S,,0^,)-T\ (28) 

Similar results for Y^'^ and Y^'^ were obtained by Fo- 
gelstrom.^ Combined with these results for the coherence 
functions, the boundary conditions for distribution func- 
tions given in Eqs. (8) and (19) provide a complete set of 
quasiclassical boundary conditions applicable to a wide 
range of non-equilibrium conditions for superconductors 
in contact with spin-active interfaces. These boundary 
conditions (Eqs. (6)-(25)) reduce to the results of Ref. 
27 for non-spin-active scattering, i.e. when Sij and S_^j 
are spin-independent.^^ 

In deriving Eqs. (6)-(25) we assumed that the inverses 
oi (3^^^ , and their partners, are defined. Equations (6)- 
(25) cannot be applied when one or more of the S matrix 
elements is zero. However, in cases where this happens 
the boundary conditions are significantly simplified, and 
can be readily derived following the procedure outlined 
in the appendix. For example, in the case of an im- 
penetrable wall, we have perfect reflection described by 
Si2 ~ S21 = 0. Then Eqs. (6)-(8) are replaced by the 
simpler set of boundary conditions. 
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rf = SnirSl, (29) 
= S,,^tsl (30) 
= Siix'^Sl^. (31) 

Similarly for perfect transmission, Sjj = 0, the boundary 
conditions simplify to 

rf = Si27M2 (32) 

ff = ^1272^^12 (33) 

= Sr2X^Sl,. (34) 

Further simplification occurs for stationary non- 
equilibrium transport. The time convolution products 
reduce to matrix products, and there are additional sym- 
metry relations: pf^ = (A^)t, pf^ = {pfjY, = (rf )t, 
i^, = (tf)t,and &tr = i<V- 

III. THE S MATRIX 

A microscopic calculation of the normal-state S ma- 
trix for a spin active interface would require a solution 
of the many-body problem in the presence of the inter- 
face potential. This is a formidable problem and outside 
the realm of a practical theory aimed at understanding 
the transport properties of heterogeneous superconduct- 
ing junctions. The alternative approach is to identify 
the structure of the S matrix, including the constraints 
of symmetry, and then model the interface in terms of 
the key physical parameters defining these characteris- 
tics, e.g. the transmission and refiection probabilities for 
normal-state electrons and holes moving along specific 
trajectories and in particular spin states. For a relatively 
small set of physical parameters, the key characteristics 
of the interface can be obtained from measurements, e.g. 
from normal-state transport properties, and then used to 
make predictions for non-equilibrium properties in the 
superconducting state. This is the most tractable ap- 
proach to interpreting and predicting the nonequilibrium 
properties of heterogeneous superconducting junctions. 

In this section we discuss the parametrization of the 
S matrix in terms of a spin-mixing angle and spin- 
dependent normal-state transmission coefficients. Other 
authors have also discussed the form of this interface 
S matrix for particular magnetic interfaces, c.f. Rcfs. 
9,10,14,32 We discuss the form of the S matrix for both 
a ferromagnetic insulating interface and a clean ferro- 
magnetic metallic interface. For both cases we assume 
the interface is atomically smooth so that the momentum 
parallel to the interface, p||, is a good quantum number. 

First consider the S matrix of a ferromagnetic insulat- 
ing or semiconducting (FI) interface. ^'"'^^''^^ Choose the di- 
rection of the spontaneous magnetization fl as the quan- 
tization axis for the conduction electron spin. Then spin 
up {+) and down (— ) electrons see the FI interfaces as 



a potential barrier with thickness I and height Eg ^ h, 
where Eg is the average band gap and h is the exchange 
energy. The effects of the FI layer on the transport of 
electrons are two fold: 1) spin filtering in which the re- 
flection (transmission) probabilities for spin up and spin 
down electrons are different, because these electrons with 
different spin polarization see different potential barriers, 
and 2) spin mixing in which spin up and down electrons 
acquire different phase shifts upon reflection (transmis- 
sion). This is the analog of circular birefringence in op- 
tics. Thus, in general the polarization of an incident elec- 
tron undergoes a rotation analogous to optical Faraday 
rotation. 

The degrees of spin filtering and spin mixing are de- 
termined by the modulus and the phase of the reflection 
(transmission) amplitudes, respectively. Consider the re- 
flection amplitude for example. In the spinor basis |±) 
which diagonalizes /t • a, the spin matrix 6*11 is diagonal. 

For an arbitrary basis 5*11 can be parameterized as 

= e*'^"/2[sii + s[,{fi ■ (T)]e*(A-)''ii/2 ^ (gg) 

where the overall phase factor, Lpu = 6^ + 9"^, and the 
spin mixing angle, = 9+ — 9-^ are defined as the 
sum and difference of the phases for the reflected spin 
up and spin down electrons. The two real amplitudes, 
sii = (|r+|-h|r_|)/2 and = (|r+| - |r_|)/2, determine 
the spin filtering effect. A similar parametrization can be 
carried out for each element, Sij, of the S matrix. 

The unitary condition, SS^ = 1, combined with sym- 
metries of the interface, provide constraints between the 
values of {(y9y , Sij, s^^}. For a specular FI interface 
with inversion symmetry, the constraint of time reversal 
symmetry, which includes the reversal of the ferromag- 
netic moment, gives ip2i = '^11 + 7r/2, and ')?2i = "i^ii- In 
this case the spin mixing angle for reflection and trans- 
mission are the same. The resulting S matrix simplifies, 
and is conveniently expressed in the basis, |±), 

c c .^2 /Rfe'''/^ \ 

52. = 5,2^.e-/^(^;"'' ^e-/^)' (37) 

where Ra + Da = 1, a =t, [■ The overall phase factor, 
ip, drops out of all observables in the quasiclassical ap- 
proximation and can be omitted. Therefore the S matrix 
is described by three parameters: the transparencies for 
spin up and spin down electrons, and and the 
spin mixing angle, 

If we also have reflection symmetry in a plane perpen- 
dicular to the interface, then S'y(— p||) = S'y (p||). This 
implies that that the S matrix for hole scattering is sim- 
ply ^ij = Sij. This model of a ferromagnetic interface 
defined by Eq. (37), as well as special cases without 
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FIG. 2: The spin up (-D+) and spin down (-D-) transmission 
probabilities, and the spin mixing angle as a function of the 
angle of incidence, $ = arccos(p/ • z). The model parameters 
for the S-FI-S interface are described in the text. 



component of the momentum of the spin down electrons 
in the FM, = [Ef cos^ $ — h)2m*, vanishes and 
becomes imaginary for larger angles of incidence. Thus, 
these spin down electrons can only tunnel through the 
FM barrier. Since the charge and heat currents are dom- 
inated by trajectories close to normal incidence we con- 
sider the transmission probabilities and spin-mixing an- 
gle in the small angle limit near normal incidence, where 
both p+ and pj are real. If we further assume the ex- 
change field is relatively weak, h Ef, then to the lead- 
ing order in h/Ef, 



9^ 
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{pt 



-p-)l (38) 
p,M)2sin2(p^0/4. (39) 



Thus, for normal incidence the spin mixing angle is of 
order d ^ pfl{h/ Ef), which can easily approach i? — *■ tt. 
The electrons of both spin species are nearly perfectly 
transmitted, ~ 1 — (/i/2£^/)^, and the spin filtering 
effect is negligibly small. 



spin filtering, have been discussed previously by several 
authors. ^'^"'"^^ The reflection and transmission probabili- 
ties, {Ra, Da, '9} are functions of the direction of the tra- 
jectory of an incident quasiparticle, p/, and depend on 
material parameters such as the band gap. Eg, the Fermi 
velocities of the electrons in the two metallic leads, -Vfi, 
exchange field, h, barrier thickness, I, etc. 

To illustrate the typical parameters for spin mixing 
and spin filtering by a FI interface consider a FI barrier 
with a band gap of Eg = 0.825 eV and and exchange 
splitting of /i = 0.18 eV, between two metallic leads. "'^ 
For conduction electrons with effective mass m* equal to 
the band mass of carriers in the FI we can calculate the 
spin-mixing angle and the transmission probabilities for 
spin up and spin down conduction electrons at normal in- 
cidence. A barrier of width / = 0.5 nm gives = 0.013, 
r>X = 0.007, and d = 0.032 tt. The ratio DJD^ vanishes 
exponentially as the barrier thickness increases, and the 
spin mixing angle -d saturates at 0.03487r, the spin mixing 
angle for a perfectly reflecting FI surface. For angles away 
from normal incidence the effective barrier thickness in- 
creases and the corresponding transmission probabilities 
decrease rapidly away from normal incidence as shown 
in Fig. 2. The spin mixing angle also decreases with the 
angle of incidence, and vanishes for grazing incidence. 

The iS matrix model in Eq. (37) is sufficiently general 
to account for the essential features of spin-active scat- 
tering by a clean, ferromagnetic metallic (FM) layer. 
For example, assume the transmission and refiection of 
electrons by the interface is controlled by Fermi wave- 
vector mismatch at the S-FM interface. Upon enter- 
ing the FM layer, the Fermi momenta for majority (spin 
up) and minority (spin down) electrons changes to p^ = 

•y/ {Ef ± h)2m* , respectively. As a result the transmitted 
majority- and minority-spin electrons acquire a relative 
phase. 

For sufficiently large angles of incidence the normal 



O 



2E 



(40) 



Thus, the dominant effect of the FM interface is spin 
mixing, and an approximate form for the S matrix of an 
ideal FM interface is given by S*!! = S22 — and 



5*12 — <5'2 



giiPgj:(ff-/i))?/2 



(41) 



Starting from Eq. (41) a more detailed model for the S 
matrix of a FM interface can be constructed by adding 
a thin nonmagnetic insulating layer, with transparency 
D, inside the ideal FM, which may describe an interfa- 
cial dielectric barrier. The composite S matrix of this 
structure takes the form of Eq. (37) with = = D. 

However, there presumably exist a wide variety of spin 
active interfaces, described by any physically allowed 
value of 1) and D^- Thus, the calculations that follow 
are carried out for a broad range of values of and Da- 



IV. FM & FI POINT CONTACTS 

We now illustrate the application of the boundary 
conditions by calculating some representative transport 
properties for both NFS and SFS point contacts. These 
calculations highlight the role of the spin mixing angle 
in modifying the local spectrum near the point contact 
and in modifying the effective transmission coefficient for 
excitations that carry currents across the interface of the 
point contact. 

Although the formalism is applicable to superconduc- 
tors with any pairing symmetry, the calculations de- 
scribed here are for spin-singlet, s-wave superconduc- 
tors. For a point contact the radius of the contact is 
much smaller than the coherence length. In this limit 
the pairbreaking effect of the FM on the magnitude of 
the order parameter can be neglected, and the voltage 
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drop occurs at the contact because of the large Sharvin 
resistance. '^^ Then at the point contact the Riccati am- 
phtudes 7^^"^ , 7^*'^^ , a; , xj^ take their local, bulk equi- 
librium values given by,^^ 



7f < A,) = -(^a2)e-'*^+^^^ 
7f > A,) = -(*a2)sgn(e)e-*^+'^ 

xf = (l-|7fntanhi^ 



(42) 
(43) 

(44) 



where we introduced the dimcnsionless parameters, 



6j = arccos — , |e| < Aj 

LXrj 



coshSj = \e\> Aj 



(45) 
(46) 



and Aj is the gap, Vj is the potential, Tj is the tem- 
perature, and is the phase of superconductor on side 
j = 1,2. Application of the boundary conditions, Eqs. 

(6)-(25), is straight forward; wc obtain T^^^, f^^"^, 

and from which wc construct the quasiclassi- 
cal Green's functions. 



A. NFS conductance 

Consider the electrical conductance of an NFS contact 
at fixed temperature, T, and voltage bias, V. Due to both 
spin-mixing and the proximity effect, the local density 
of states (DOS) of the superconductor deviates from the 
bulk BCS form. Surface states appear below the gap, and 
a splitting of the DOS for spin up (+) and spin down (— ) 
excitations develops for any < < tt. 



iV±(|e| < A) = 7V;- 



1 - RfRi 



N±i\e\>A)^Nf 



R^Rl - 2y/R^icos{2S T ^) 

g25 ^ e-^^R^Ri - 



(47) 



(48) 



where Nf is the density of states at the Fermi level, and 
6 is defined in Eq. (45) and Eq. (46). For perfect reflec- 
tion, i?! = i?| = 1, there is a true surface bound state, 
analogous to the Shiba state bound to a magnetic impu- 
rity in an s-wavc superconductor.'^^ For finite transmis- 
sion, the surface bound states broaden into resonances 
due to the proximity coupling with the normal metal. 
In the tunnelling limit, i.e. for low transmission with 
Z)| ~ Z3| = 13 ^ 1, N±{€) exhibits a relatively sharp 
resonance peak below the gap at e± ~ ±Acos(i9/2) with 
a width of order 7 ~ DA/2.^^ For higher values of the 
transmission probability the resonances broaden into a 
sub-gap continuum. 

The differential conductance for low-transmission junc- 
tions reflects the resonance states which transport charge 
via resonant Andrcev reflection. The spectral current 
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FIG. 3: Differential conductance, dl/dV, for NFS contacts 
with different spin mixing angles, 1? = 0, 7r/2, vr at T = and 
for Z?! = D| = D = 0.05 and D = 0.1. The normal-state 
resistance of the point contact is Rn = {e^ NfD A)~^ , where 
A is the area of the contact. 



density, j(e), can be calculated from the solution for the 
quasiclassical propagators at the interface. 



3N 

< 



AcoshS[e-^D^Di/{Df + D^) + sinh(5] 



'- = 2y- 

JN ^ 1 



R^Rl - 2y%R7cos(2(5 ± t9) ' 



(49) 
(50) 



where j< is the spectral current density for |e| ^ A and 
Jn oc e{D^ + Dl) is the spectral current density when 
the both electrodes are in the normal state. The total 
current density is then given by 



] 



\ J rfej(e)[tanh 



eV 



2T 



tanh 



eV 



2T 



(51) 



Figure 3 shows the zero temperature differential conduc- 
tance for NFS junctions with different spin mixing an- 
gles. The proximity effect is evident as a finite sub-gap 
conductance even for a non-magnetic interface. The in- 
terface resonance induced by a finite spin-mixing angle is 
also clearly exhibited as a broad peak in sub-gap conduc- 
tance at eV « Acos(i?/2). Note also that the width of 
the resonance, 7 « DA/2, provides a spectroscopic mea- 
sure of the interface transmission probability. However, 
thermal broadening dominates the width of the sub-gap 
resonances in the tunneling limit, except at very low tem- 
peratures, as shown in Fig. 4. 

Asymmetry in the transmission probabilities, ^ 
D)^, for spin up and spin down excitations also leads to 
a finite spin current. The corresponding spectral current 
density is given by 



2(L>I - Di) sinh(2(5) 



i?Ti?X~ 2 /Rp?r COS 



(52) 
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FIG. 4: Differential cliarge conductance, dl/dV, and spin 
conductance, dP /dV, for an NFS contact witli spin mixing 
angle, ■& — 7r/2 and transparencies, Z)| = 2_Dj — 0.2, at 
temperatures, T = O.OlTc, T = O.lTc and T = 0.2rc. Note 
that the voltage is normalized in units of Ao, the gap at T = 0. 



for |e| > A. In contrast to the charge current the sub-gap 
spin current spectral density vanishes identically (j*^ = 
0) because there is no resonant Andreev reflection for spin 
transport. The total spin current is then given by Eq. 51 
with j{e) — * i^le). In Fig. 4 we also show the differential 
spin conductance for weak spin filtering, Di = 2Di = 
0.2. Note the onset of the spin conductance at eV = A 
for T — > 0, and the absence of Andreev resonance peaks 
in the spin conductance for non-zero spin mixing. The 
spin conductance is normalized by the normal-state spin 
conductance of a point contact of area A, (i?^)^^ = 
NfiD^-D^)A. 

The limit of extreme spin-filtering provided by a half- 
metallic ferromagnetic metal is discussed in detail using 
the transfer matrix method to incorporate spin-mixing 
at the interface between a FM and a superconductor.^"' 
We obtain the same results as Ref. 24 for the current of 
a half-metallic ferromagnetic-superconductor point con- 
tact by setting = D and — 0. In this limit the 
sub-gap conductance from resonant Andreev reflection is 
completely suppressed by the spin-filtering effect. 



B. SFS thermal conductance 

In a recent report we discussed the role of Andreev 
bound states in regulating quasiparticlc transport of heat 
through point-contact Joscphson weak links. '^^■^^ Spin 
mixing generates a spin-resolved spectrum of Andreev 
bound states at a spin-active point contact even in the 
absence of a phase bias. We apply the formalism and 
boundary conditions developed in previous sections to in- 
vestigate the effect of spin mixing on the phase sensitive 
heat transport in temperature biased SFS point contacts. 
We show that the relative phase shift of spin up and down 
electrons, together with the phase bias 4> = 4>2 — 4>i, deter- 



mines the spectrum of Andreev bound states at the point 
contact. The effects of these states on the transmission 
probability of continuum excitations that transport heat 
is calculated. 

The thermal conductance of the point contact is de- 
fined by the ratio of the total heat current and the tem- 
perature bias 6T = T2 — T1 in the limit 5T — !■ 0. The re- 
sults reported below are normalized by the normal-state 
thermal conductance at Tc- — j^ANfVfTc{D'^ +Di), 
where A is the area of the point contact. 

Following the same line of argument as described for 
the NFS contact we apply the boundary conditions, Eqs. 
(6)-(25), to construct the Green functions for the SFS 
contact. The ABS spectrum is straightforward to calcu- 
late and has been discussed in the context of the Joseph- 
son current of the SFS weak link by several authors. ^'""^'^^ 
There are two branches (labelled as ±) of spin-up An- 
dreev bound states with energies. 



ei^ = A sgn sin 



'd±p 



where the angle p is defined as 

p — arccos( ^/R^ Ri + ^ cos i 
The spin down bound states are at energies 



(53) 



(54) 



(55) 



At i5 = or TT, the ABS spectra is degenerate with respect 
to spin. For < ■)? < tt the spin degeneracy is lifted, thus, 
the typical bound state spectrum has four branches, two 
branches per spin direction for p ^ 0. Branches with 
opposite spin are "mirror reflections" of one another with 
respect to the Fermi energy (e = 0). In the tunnelling 
limit, D| = D| = D — > 0, and to leading order in p ^ t?, 
the splitting of the spin up states is given by 



1 ~ A[cos^Tsin^%/:Dsin^]. (56) 



The spectral weight of an ABS comes at the expense of 
the continuum spectrum (|e| > A). In addition there is 
asymmetry with respect to spin, A^^(|e| > A) 7^ iV-l-del > 
A), and as a result spin up and spin down quasiparticles 
contribute differently to the heat current. To compute 
the heat current, we follow the procedure described in 
Ref. 33. The boundary conditions for the distribution 
functions, Eqs. (8) and (19), enable us to obtain an an- 
alytical result for the Keldysh Green's function at the 
point contact. From this result we can calculate the spec- 
tral density for the heat current for the set of trajectories, 
{li, lo, 2i, 2o}. The spectral heat current contains con- 
tributions from spin up and spin down electron-like and 
hole-like quasiparticles. 



p/) ^ ^Tr[5^^ - gf„] = ^Tr[5f„ 



9l\ 



(57) 



where Tr represents the trace of Nambu matrix propaga- 
tors. It is straightforward to show that the spectral heat 
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current can be expressed in an intuitive form by intro- 
ducing the effective transmission coefficient S) for heat 
transport, 



/(£,P/) = 



-2e tanh 

2Ti 



e 

- tanh 

2T2 



Vie,pf). (58) 



The contributions to X) come from direct (e — > e) trans- 
mission and branch conversion {e{h) hie)) transmis- 
sion channels, D = Se^e + 1)e^h; 



+L»j.)cosh(2(5) - 2^Di^Di cos i? cos 0)] 
xsinh^V^^ (59) 



[{D^ + i?i) - 2D^D^ - 2^R^RiD^Di cos0] 
xsinh^V-Z^ (60) 
[v/i?T^i + V-^T^i cos <f) - cosh(2(5) cos 1}] 2 



+ sinh^(25) sin^??. 



(61) 



In the normal state Xl [D^ + Z3|)/2, and spin mixing 
has no effect on the quasiparticlc transport. However, 
in the superconducting state the transmission coefhcient 
D for quasiparticles of energy e and momentum p/ is 
sensitive to both the phase bias, (p, and the spin mixing 
angle, §. 

Consider first the case with (p — Q. Spin mixing leads 
to bound states at ee = ±A cos('!9/2). For the case in 
which D-t = ~ D, only the direct transmission chan- 



nel contributes, i.e. De 



= 0) = and 



0) = D- 



A2 C0S2 I 



<D. 



(62) 



Thus, quasiparticle transmission is suppressed by the 
spin mixing effect for any value of the normal-state trans- 
parency, D, and for all energies. The suppression is most 
severe a.t -d ~ n, i.e. when the ABS is more strongly 
bound. ^'^ Thus for (j) = spin mixing suppresses heat 
transport for any value of the normal state transparency. 

For a spin-inactive point contact (?? = 0) in tunnelling 
limit, D ^ 1, it is known that 5)(e) has a resonance peak 
at 

eres/A ^l + ^D sin2 ^ + 0{D^) , (63) 

which is a reflection of a shallow bound state just below 
the gap edge at eb/A = f - \D sin^ | + 0{D^). Tuning 
(j) from to TT leads to an increase in the thermal conduc- 
tance because of resonant transmission of quasiparticles 
Enhanced transmission still exists for SFS 



at e 



33 



point contacts, but as d increases the resonance peak of 
S)(e) gradually vanishes. For i9 = tt, the bound states 
are at energies eb = ±A-\/Z? sin(0/2), and there is no 
resonance peak of Ti{e). Instead D is suppressed at all 
energies. 



D- 



2- AZDsin^l 



<D. 



(64) 



Thus, to leading order in D, D{-d = tt) ~ D(l - A2/e2) 
is independent of 0, so phase modulation of the thermal 
conductance vanishes. These features are shown clearly 
in Fig. 5 for the thermal conductance calculated in the 
tunnelling limit with = = D = 0.1 and for gen- 
eral values of and (j). Resonant enhancement of the 
conductance occurs in the vicinity of sa tt and 1} 1. 
Increasing suppresses the overall thermal conductance 
as well as the phase modulation. 

Figure 6 shows the thermal conductance in the high 
transparency limit with = = D = 0.9. At = 
(a "0" junction), tuning towards tt pushes the ABS 
deep into the gap, so D{e) is increasingly suppressed from 
D, and the thermal conductance goes down. The oppo- 
site occurs at = TT (a "tt" junction^): tuning (f> from 
to TT pushes the ABS toward the gap edge, so the ther- 
mal conductance increases. The phase modulation of the 
thermal conductance for a general value of d can be un- 
derstood qualitatively in a similar manner. The thermal 
conductance is maximum when the bound states are clos- 
est to the gap edge. The different phase modulation of 
the thermal conductance for "0" junctions (i.e. < tt/2) 
versus "tt" junctions (i? > tt/2) should be observable in 
high transparency SFS junctions; one should in princi- 
ple be able to change the spin mixing angle by varying 
the thickness of the FM layer, thus tuning between "0" 
or "tt" junction behavior. The phase of the SFS junc- 
tion can then be controlled by varying the magnetic flux 
linking a SQUID containing the SFS contact. 

In contrast to FM contacts, SFS junctions with FI 
contacts are expected to be in the tunnelling limit, i.e. 
D-^^l <^ 1. For the FI interface described in section III 
with / = 0.5 nm, D| = 0.013, = 0.007 and d = 
0.0327r, the discrimination between transparencies for dif- 
ferent spin orientation is relatively large, ~ 20^, but 
the spin mixing is weak. As a result the phase modulation 
of the thermal conductance is almost the same as that of 
a spin-inactive tunnel junction. '^^ Figure 7 shows a map 
of the thermal conductance as a function of temperature 
and phase bias for the junction parameters given above. 
For simplicity wc neglected the angular dependence of 
Da and d. Note that resonant transmission leads to en- 
hancement of the thermal conductance compared to the 
normal-state conductance at Tc over a broad temperature 
range 0.5Tc < T < for (/) > n/2. 



V. DISCUSSION AND CONCLUSION 

In this paper we derived the boundary conditions for 
the quasiclassical Ricatti amplitudes, including the dis- 
tribution functions, , for general interfaces. This com- 
pletes the development of the boundary conditions for 
the Riccati amplitudes at spin-active interfaces initiated 
in Refs. 9. The results summarized in section II apply 
to a broad range of superconducting interfaces, and are 
specifically applicable to non-equilibrium transport in- 
volving spin-active interfaces. The boundary conditions 
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FIG. 6: The thermal conductance for an SFS point contact 
with Dt ^ Di= 0.9 at T = O.STc. 



can be applied to investigate dynamical properties, are 
applicable for any paring symmetry. The effects of disor- 
der are readily described within the quasiclassical theory 
framework. 

In the dirty limit the quasiclassical theory for 
conventional superconductors can be formulated in 
terms of Fermi surface averaged quasiclassical Green 
functions. '^^ The quasiclassical transport equations re- 
duce to diffusion-type (Usadel) equations. Boundary 
conditions for the the Fermi-surface averaged Usadel 
propagators at a spin active interface has been discussed 
recently in Ref. 38. A comparison between results for 
the dirty SFS and NFS junctions based on the Usadel 



equations and boundary conditions and the dirty limit of 
the Ricatti formalism with the general spin-active bound- 
ary conditions developed here will be discussed in a sep- 
arate report. Finally, we note that we have presented 
boundary conditions and results applicable to specular 
interfaces. The boundary conditions can be extended to 
include atomic scale roughness at a spin-active interface, 
but this extension is outside the scope of this report. 

The Riccati approach to the boundary problem in 
quasiclassical theory is complimentary to the T-matrix 
approach, ^^'^-'^"^^ which has been used to calculate the 
sub-harmonic structure of the current-voltage charac- 
teristics of magnetic Josephson point contacts, and 
recently to study the Josephson effect and conduc- 
tance of superconductors in contact with half-metallic 
ferromagnets.^^'^* The advantage of the Ricatti approach 
is that the implementation of the boundary condition is 
given explicitly in one step, without introducing auxiliary 
propagators. The Riccati formulation also makes analyti- 
cal analysis more tractable. The Riccati method has been 
successfully used to study the equilibrium properties, e.g. 
the d.c. Josephson effect and temperature-induced — tt 
transition of SFS junctions. ^'^^ As shown by Barash et 
al,^^ the Riccati approach is especially powerful to study 
more complicated hybrid structures such as SFIFS junc- 
tions. 

As examples of the potential application of the newly 
derived boundary condition for the non-equilibrium dis- 
tribution functions, , we investigated the charge and 
spin conductances for NFS point contacts and heat trans- 
port through temperature and phase biased SFS point 
contacts. We showed how the spin-mixing angle iD, de- 
fined as the relative phase between spin up and down 
electrons upon transmission (or reflection), controls the 
local excitation spectrum and the transport of charge, 
spin, and energy across the point contact. Beyond this 
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relatively simple application, we expect that new results 
and novel physics can be explored for a broad range of 
non-equilibrium transport problems involving spin-active 
interfaces with the new boundary conditions. 
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APPENDIX: DERIVATION OF EQS. (6)-(25) 

In this appendix we outline how the boundary condi- 
tions for the Riccati amplitudes, Eq. (6)- (25), are derived 
from the Millis-Rainer-Sauls boundary conditions. Our 
starting point is the set of Eqs. (63)-(66) of Ref. 26, 
which in terms of Shclankov projection operators, '^^ 

(A.l) 

can be written as 

Pl''iS22Pl'Sl^ + S2iP^'Sl,)Pl° = , (A.2) 



P+{Sl2P^°S22 + S\^P+Si2)P^- = , (A.3) 
Pl\S\^Pl°Sii + Sl^Pl°S2i)P}' = , (A.4) 
+ Si2P+S\^)P+ = . (A.5) 

Following the convention in literature, we denote a 
Keldysh(8)Nambu matrix with check accent, a Nambu 
(particle- holc(3)spin) matrix with hat accent, and a spin 
matrix without any accent. We use the ansatz of 
Eschrig^^ to parameterize the projectors in terms of Ric- 
cati amplitudes, e.g. for we have 

P^ - ( J^)(l-7fffri(l, 7f) (A.6) 
Pt = ['\'){l-lt^ir\lt. 1) (A.7) 

p^ - p^{i (o if )^+^^-^) 

where we omitted the superscript \i. The strategy to 
solve Eqs. (A.2)-(A.5) is to reduce the order of the equa- 
tion set by exploiting the properties of the projectors and 
to decompose the equations for the Keldysh(g)Nambu ma- 
trices into equations of spin matrices. 

Take the retarded components of Eqs. (A.2)-(A.5) and 
plug in the expressions for P_|_' . Each of the four equa- 
tions of Nambu matrices collapses into an equation for 
the coherence functions which are spin matrices. For ex- 
ample, the retarded part of Eq. (A.3) and (A.4) give 



- rf 7f )-'(rf - 51272'') - {sl^^^l - 12^2-^iX - 72''rf )-^/32«2, (A.9) 

/^fi(l - rf 7f )-^(rf - 5n7f ) = {S\^n - 7f 5ti)(l - ~^^Y^r^~p^2- (A.IO) 
Eqs. (A.9) and (A.IO) suggest the solution 

(rf - 5n7f )(/3f2)-' = (rf 5i2 - 5i272^)(/32^2)-\ (A.ll) 

(/3fi)-i(5tirf - 7f^t J = (/3f2)-'(5l2r^ - 1^3^22)- (A.12) 



The above results assume the inverses of (5^ and their 
~ partners exist. This is generally the case for partially 
transmitting interfaces. Otherwise, Eqs. (A.9) and 
(A.IO) can be solved rather trivially since some of the 
S matrix elements vanish. Equations (A.11)-(A.12) im- 
mediately lead to the result obtained by Fogelstrom,^ Eq. 
(26), which is equivalent to Eq. (6), as well as 

= r2l72''^L+i^7f^k (A. 13) 

= +Pf2)~'52^2-(/?fi)-'52^i]-'(/?f2)-'(A.14) 

tl = -[{0^2r'si2'iP?ir'si,]-\p^,r\A.i5) 

It is also straightforward to verify that the solution, Eq. 



(26) and Eqs. (A.13-A.15), indeed satisfy Eqs. (A.9) and 
(A.IO). In a similar manner, the boundary conditions 
for all retarded and advanced coherence functions can be 
derived. 

The Keldysh components of Eqs. (A.2)-(A.5) can be 
simplified by using the equations for retarded and ad- 
vanced projectors. Then after plugging in the expressions 
for all the projectors, p^^/^/^ ^ once again we find that 
the Nambu matrix equations collapse into spin matrix 
equations for the Riccati amplitudes. For example, the 
Keldysh components of Eqs. (A.3) and (A.4) lead to 
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+/3fi(l - rf7f (1 - 7i^f^^)-iai\ 

-/32« (1 - rf7f )-i (1 - 72^f^)-'«^i 

-{Sir? - 7f^L)(i ~ 7f rf )-i if (1 - ftjfr\rts^, - ^1,71^) 

+ (^^ir? - 7f^L)(l - 72^rf )-i if (1 - f^72^)-Hf ^521 - ^2i7i^), (A.16) 

= [-/?f2(i - rf7f )-i^i2 + isl,r^- - 72''-^L)(i - 72''rf )-i72''^22] 
+^2^ [-^l27i^(i - th^rHrts,2 - ^1272^) + ^2^2(1 - 72^r^)-'«^2] 

+Pl2(l-il7l) -^1 (l-7lilj "12 

(1 - r^^'^)-' ^2^ (1 - 72^f ^)"'«4 

-{SUr?^ - 72^^^l2)(l - 7f rf )-i if (1 - ft^f)-\ftSu - ^1272^) 

+ (^l2rf - 72''^^L)(1 - "fMr' if (1 - f ^72^)-'(f ^^22 - ^2272'')- (A.17) 



Again, if the inverses of P^j do not exist Eqs. (A.16) 
and (A.17) simpUfy and are readily solved. For the gen- 
eral case, in order to solve for and , we construct 
(/?2«i)-i X (£;g. A16) X (afi)-i - (/Jg)"' x {Eq.A.lT) x 
(ai,)-' and (/Jfi)"^ x (i?^. A.16) x (a^\)-i - (/Jf^)"! x 
(i?^. A.17) X {a^2) ^ to obtain transformed equations, 



Eq. (A. 10)' and Eq. (A. 11)', which are not repro- 
duced here. The transformed Eq. (A. 10)' contains only 
,X2 and X2 ■ We regularize each term by car- 
rying out a series of transformations. For example, the 
terms proportional to xf are transformed to 



+ (/32«i)-^/3fi(l " rf 7f )-^(af 72^5.1 - Sn) {a^)-' + (P^,)-' xf (5,^72^ - - 7i^f ^^)-^<K)-^ 

+ iP^i)-' xUaAr\ (A.18) 
and further into 



(/32^i)-'/3fi(l - r?'J^r'r?,{Sl)-'x^{Sn)-'rUl - j^fty'atAa 



+ iP^2)-'f3?2{l - r?l^r'rf,{Sli)-'MSii)-'rt{l - l^rt)''at2ia^2r' , (A.19) 
with the help of algebraic identities such as 

a?i l^S2i +r^i = Sn, a^i ^^822 + tf, = 8,2, (A.20) 
l-rf7f = ^H/5fi+iM, (A.21) 
rf 7f )-4f /?f2 = {m-'P?2{^ - rf 7f )-V^^(/3fj . (A.22) 

As a result the transformed Eqs. (A. 10)' and (A. 11)' can be expressed as 

(/3« )-^/3fi(i - rf7f )-ini(i - it^^r'PtM,)-' = {P^2)-'PUi - rf7f )-ini(i - ^t^tT' pUp^2)-\{^-'^^) 

{P?2)-'PU^ - T^l^r'n2{l - 1^^)-'PUP21)-' = (/3fi)-'/32'i(l - rf72'')-'"2(l - 72^f^)-'/?2l(/3n)-\(A.24) 



with or cquivalcntly, the original Eqs. (A.16) and (A.17). The 

equation ni = yields the boundary condition for Xf 
in Eq. (8). The boundary condition for Xf is derived in 



n, ^ r^.x^Pi - tf^x^tt + a^ii^c^t , (A.25) 



712 = X^ T2iX2^2r ~ ^2;^f ^2r + 0'2iXi o.2r ■ (A. 26) a similar manner starting from Eqs. (A. 2) and (A. 5). 



Obviously m = 712 = satisfies Eq. (A. 23) and (A. 24), 
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